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Effective field theory is considered to provide a highly useful framework for connecting 
nuclear physics with the symmetries and dynamics of the underlying theory of strong 

t— i ■ interactions, QCD. Of many issues that are of great current interest in this domain, I 
concentrate here on two: (1) A new class of ab initio calculations of observables in two- 

i/-) \ nucleon systems; (2) Attempts to extend chiral perturbation calculations to higher-order 
terms. 

o : 

0\ ■ 1. Introduction 

5 : 

' One of the major challenges in nuclear physics today is to establish a connection between 

nuclear dynamics and the fundamental QCD. Effective field theory (EFT) provides a 
natural and useful framework for this purpose. The basic idea of EFT is simple [0,0. 
Suppose we are interested in phenomena characterized by a typical energy-momentum 
scale Q. We expect that the degrees of freedom whose energy scales are significantly 
larger than Q need not feature explicitly in our Lagrangian. So, introducing a cut-off 
scale A that is reasonably large as compared with Q, we separate our fields (generically 
denoted by $) into a high-energy part $ H and a low-energy part $ L . Integrating out <E> H , 
we arrive at an effective Lagrangian that only involves $ L . The original Lagrangian £ 
and the effective Lagrangian C eS are related as 

J[d$]jS* xC W = [[d^e'S^^ (1) 

By construction, £ e ff inherits all the symmetries (and the patterns of symmetry breaking, 
if any) of the original C. Then £ e g- is given as the sum of all possible monomials of $ L 
and their derivatives that are consistent with the symmetry requirements. Since a term 
involving n derivatives scales like (Q/A) n , we have perturbative expansion with respect 
to d^/A. In a nuclear-physics application of EFT, the original Lagrangian £ is the QCD 
Lagrangian, but, if we are concerned with the energy- momentum regime Q <C A QCD ~ 1 
GeV, the relevant effective degrees of freedoms are not quarks and gluons but hadrons. 
Furthermore, for Q < m n , it is reasonable to retain only the pions and nucleons as 
effective degrees of freedom. Meanwhile, chiral symmetry of QCD must be respected in 
our effective world. The resulting effective theory is called chiral perturbation theory 
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(%PT) fHH. In fact, the inclusion of the nucleon in %PT poses a problem because its 
mass is comparable to A QCD . The heavy-baryon chiral perturbation formalism (HB^PT) 
allows us to circumvent this difficulty £™ has, as effective degrees of freedom only 
the pions and the large components of nucleons, and it involves expansion in <9 M /A QCD , 
m 7r/A QC D and d^/m N . Since m N ~ A QCD , we usually lump together chiral and heavy- 
baryon expansions. In this combined expansion scheme, the effective chiral Lagrangian 
can be organized in terms of the chiral order index v defined by v = d + (n/2) — 2, where 
n is the number of fermion lines that participate in a vertex, and d is the number of 
derivatives (with m n counted as one derivative). The leading order terms are given as || 

f 2 

£ (o) = iJL^ld^d^U + ml(U ] + 17-2)] 

+B(iv .D + g A S-u)B-W C A (BT A B) 2 (2) 

1 A 

£(i) = -^B{S-D,vu}B + 2 Cl mlBBTi(U + -2) 
2m N 

+ (c 2 --^)B(vu) 2 B + c 3 Bu-uB 
8m N 

( ' 9 (BB)(BiS-uB) - —(BS^B^Bm^B). (3) 



2m N 2m N 

Here U(x) is an SU(2) matrix field related non-linearly to the pion field, £ = y/U, = 
i(^®iJ>£ ~ £>n = ij5< J Liv vU an d 77 M is the covariant derivative acting on the 

nucleon. Here we have retained only those terms that are relevant to our subsequent 
discussion. The counting rule of Weinberg || is that a Feynman diagram consisting of 
A nucleons, N E external fields, L loops and C-separated pieces is order of 0(Q V ) with 
v = 2L + 2{C - 1) + 2 - (A + N E ) + v t . 

In fact, straightforward chiral counting fails for a nucleus, because purely nucleonic 
intermediate states that occur in a nucleus can have very low excitation energies which 
invalidate the ordinary chiral counting rule ||. To avoid this difficulty, we classify Feyn- 
man diagrams into two groups. Diagrams in which every intermediate state contains at 
least one meson in flight are called irreducible diagrams, and all others are called reducible 
diagrams. The chiral counting should be applied only to irreducible diagrams. The con- 
tribution of all the irreducible diagrams (up to a specified chiral order) is then used as an 
effective operator acting on the nucleonic Hilbert space. By summing up infinite series of 
irreducible diagrams (either solving the Schrodinger equation or the Lippman-Schwinger 
equation), we take account of reducible diagrams. This two-step procedure may be re- 
ferred to as nuclear chiral perturbation theory [|3]||. 



2. Hybrid approach to nuclear x?T 

In applying nuclear xPT to cases that involve external probes, a nuclear transition 
operator T is identified with a set of all the irreducible diagrams (up to a given chiral 
order v) with an external current inserted ||. In a fully consistent xPT calculation, 
T is to be sandwiched between the initial and final nuclear states which are governed 
by the nucleon interactions corresponding to the z/-th order irreducible diagrams. In 



3 



practice, however, we often use initial and final nuclear wavefunctions obtained from the 
Schrodinger equation that involves phenomenological nucleon-nucleon interactions. This 
eclectic treatment may be called a hybrid approach to nuclear xPT. The hybrid approach 
was used extensively to test the validity of the "chiral filter mechanism" |7|||. This 
mechanism is the statement that soft-pion exchange (unless suppressed by symmetry or 
kinematics) should give a dominant exchange- current contribution. Rho || gave a clear 
interpretation of this dominance in the language of chiral counting. As argued in j7|, the 
space component of the vector current (V) and the time component of the axial current 
(A ), in general, are expected to exhibit this dominance. Regarding the V, a detailed 
calculation based on the hybrid HB%PT was carried out by Park, Min and Rho §| for 
the isovector Ml transition amplitude in the n(thermal) + p — > d + 7 reaction. The result 
of their next-to-next-to-leading order calculation indicates that the soft-pion exchange 
current indeed gives a dominant contribution and that, with the next order corrections 
added, agreement with the data is perfect. As for the A , the enhancement factors for 
first-forbidden (AJ = 0, An = yes) /5-decay transitions were calculated in the hybrid 
HB^PT [[10 1 . The results indicate that the bulk of the empirical enhancement factors |H| 



can be attributed to the soft-pion exchange current, and that the next-order corrections 
are a small fraction of it. 

Furthermore, Park, Min, Rho and myself (PKMR) |12| have recently computed in 
HB^PT the cross sections for the solar proton burning process p + p — > d + e + + u e . The 
pp fusion rate obtained in PKMR's calculation essentially agrees with the rate used in the 



standard solar model fll3|. This result is important in the light of Ivanov et al.'s recent 



claim [14] that their field theoretic approach gives a pp fusion rate significantly different 



than that obtained from the ordinary nuclear physics approach based on Schrodinger 



equations. PKMR's result based on xPT does not support the claim in Ref.[14 



3. First attempt at ah initio calculation 

Despite the impressive success of the hybrid approach, it is important to formally justify 
or go beyond this approximation. Very recently we (PKMR) have carried out a calculation 
in which the transition operators and the nuclear interactions are treated on the same 



footing ||I5|,|IEf. This formally consistent treatment may be called an ah initio calculation. 



Let me summarize here some salient features of our work described in |16 . 

We have performed next-to-leading-order (NLO) ah initio calculations for the two- 
nucleon systems both with and without the pion field. The purpose of considering these 
two cases is to examine whether the "pionless" and "pionful" effective theories, which 
have different cut-off scales, exhibit behaviors that are generally expected for EFT. 

Since, to NLO, pion loops do not enter, we can simply work with a potential. The bare 
potential V has the form 

V(q) = -n ■ r 2 A ^ ' ^ + — [C + {C 2 S» + D 2 ^W] q 2 (4) 

with a ij = 3/V8[(a{oi + a{a l 2 )/2 - {8 ij /3)a x ■ a 2 ], where q is the momentum transfer. 
The parameters C's and D 2 are defined for each isospin channel. The first (nonlocal) 
term is the pion exchange involving the Goldstone boson and hence completely known 
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from xPT. The (local) terms in the square brackets represents the effects of the degrees 
of freedom that have been integrated out. Use of the potential Eq.(|4]) in the Lippman- 
Schwinger equation generates an infinite series of divergent terms. How to regularize this 
divergence is one of the hottest issues in nuclear xPT (see blow). We use a momentum cut- 
off scheme |L7j and introduce a Gaussian cutoff: V(r) = J[<i 3 q/(27r) 3 ] e iqr ^(q 2 ) V(q), 
with S*A(q 2 ) = exp(— q 2 /2A 2 ). For a given cutoff we can determine the constants C's 
and D's in (f|) by relating them (after renormalization) to the scattering length and the 
effective range for the scattering channels or to a selected set of the deuteron observables. 
Then we are in a position to predict the N-N scattering phase shifts and the low-energy 
properties of the deuteron (other than those used as input). Moreover, we can make 
parameter-free estimation of electroweak transition amplitudes. (In our NLO calculation 
this is limited to the 1-body (impulse- approximation) contributions.) For this we only 
need the one-nucleon electroweak current derived to leading order. We have evaluated 
one-body contributions to the charge radius r d , the quadrupole moment Qd, and the 
magnetic moment \id of the deuteron. Furthermore, we have computed the one-body Ml 
matrix element M. M \ for the np capture process and the Gamow- Teller matrix element 
Ai GT for the pp fusion process. The upshot of PKMR's ab initio calculations is as follows. 
(1) All the calculated quantities are in good agreement with the empirical information in 
both pionless and pionful cases. (2) The results are stable against the variation of the 
cut-off parameter A, so long as it lies within a reasonable range; this reasonable range 
is found to be A =100 - 300 MeV (A =200 - 500 MeV) in the absence (presence) of the 
pion, values that are consistent with the general EFT consideration. (3) The presence 
of the pion brings a noticeable improvement in the accuracy of the prediction, markedly 
reducing the cutoff dependency. 

The above-mentioned "good agreement" requires a little more explanation. First, for 
the quantities that do not involve external probes, our calculated values can be compared 
directly with the experimental values, or with the values obtained from e.g. the Argonne 
v 18 potential (which was constructed to fit the data). The experimental phase shifts are 
reproduced very well up to p ~ 70 MeV without the pion and up to p ~ 140 MeV with 
the pion. The deuteron D-state probability corresponding to the vl8 potential is also 
reproduced satisfactorily. People may ask, however, what is the difference between the 
familiar effective-range formula and our prediction. It is true that, in the pionless case, 
the two low-energy constants, C and C 2 in Eq.(g), practically replace the roles of the 
effective-range expansion parameters, a and r e . With introduction of the pion, however, 
there is no such trivial correspondence, and xPT expansion contains more physics in 
it. Furthermore, xPT enables us to apply a unified expansion scheme not only to N-N 
scattering but also to nuclear transition processes. 

The quantities that involve electroweak probes may be classified into two kinds accord- 
ing to whether two-body contributions are expected to be very small or rather significant. 
To the first kind belong r d , Q d and fid, while A4 M1 and AA GT are of the second kind. 
As for the first group, PKMR's results with the one-body current alone agree well with 
the experimental values. As for M. M1 and Ai GT , PKMR's results are found to be in 
good agreement with the one-body contributions calculated with the use of v 18 poten- 
tial. This aspect is of course welcome but precise comparison with the experimental A4 M1 
can be done only after inclusion of exchange-current contributions. (For A4 GT there is 



5 



no experimental data.) Since a completely self-consistent evaluation of the exchange- 
current, which requires an NNLO calculation, is yet to be done, we may proceed as 
follows. We introduce the ratio TZ of the 2-body contribution to the 1-body contribution 
by M M1 = .A/C ody + -A^Mi° dy = A^Mi° dy (l + TZ). Park, Min and Rho @ calculated TZ in 
hybrid HB^PT with the use of the vl8 potential and the resulting A4 M1 showed excellent 
agreement with the experimental value. Since the ratio TZ is very likely to be much less 
sensitive to the wavefunctions than Ai M1 itself, it is reasonable to expect that TZ calcu- 
lated in the hybrid approach is close to TZ that would result from an ab initio evaluation. 
This expectation gets additional support from the fact that the range of the soft-pion 
exchange current is about the same as the range probed in our EFT. In this sense the 
hybrid approach is justifiable from the EFT point of view. 

In the above we introduced a momentum cut-off to regularize the divergence that ap- 
pears in the Lippman-Schwinger equation with a potential of the type Eq.(^). Another 
regularization method, called the polynomial divergence subtraction (PDS), was proposed 
by Kaplan, Savage and Wise and it has been used extensively to calculate various 
observables in the two-nucleon systems |fL9l,j20|l. One of the advantages of PDS is that it 



preserves chiral invariance. By contrast, the cut-off regularization loses manifest chiral 
invariance. Although this feature is not consequential in the NLO calculation discussed 
above, it does become relevant in higher-order calculations. Meanwhile, the PDS count- 
ing scheme generates a counter term that is of the same order as the leading-order pion- 
exchange current |20 |. Since the strength of this counter term needs to be determined 



using data, there is significant loss of predictivity. This is to be contrasted with the cut- 
off scheme in which the leading-order pion-exchange contribution is a prediction of the 
theory, and in which the problem of unknown counter terms occurs only in the loop cor- 
rections to the one-pion exchange term [|,1C]. For more discussion on the regularization 
schemes in the two-nucleon systems, see e.g. |f2T|. 



4. p + p^p + p + n reaction near threshold 




Fig. 1(a) 

I next discuss the near-threshold pp 




rO 



tracting a great deal of attention \T1 

study this very specific process. Besides the availability of high-precision data [p2 



ppn" reaction, which has recently been at- 
91. You may wonder what motivates us to 

a 



strong motivation comes from the exceptional sensitivity of this reaction to higher chiral- 
order terms. The corresponding charged-pion production reaction, e.g. pp — > pnn + , is 
described reasonably well by the single nucleon process (Born term), Fig. 1(a), and the 
"large" Weinberg- Tomozawa s-wave pion rescattering process, Fig. 1(b). By contrast, the 
Weinberg- Tomozawa term does not contribute to the pp — » ppir reaction, rendering it 
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particularly sensitive to and hence a good testing ground for the less-well-understood 
"small" isoscalar s-wave pion rescattering terms. 
The first ;^PT calculations for pp — > ppn° were carried out by Park et al. [25| and 



Cohen et al. ||26|| . In the hybrid HB^PT they used, the transition amplitude is given by 
T = ($/| T^'l&i), where (|$/)) is the initial (final) two-nucleon state distorted 
by the phenomenological N-N interaction. T^ v ' stands for the transition operator of 
chiral order v. The lowest-order one-body impulse term [Fig. 1(a)] gives T^ v= ~ x \ while 
the lowest-order two-body rescattering term [Fig. 1(b)] gives T^ u=1 \ The u — 1 HB^PT 
calculations in f25| , |2"6|| lead to pp — > pp7r° cross sections that are much smaller than the 
measured values. Meanwhile, Lee and Riska's work p3| based on the one-boson exchange 
N-N potential suggests that shorter range isoscalar meson-exchange processes, like a- 
and uj- exchanges between the two protons, might be very important for the pp — > pp7r° 
reaction. To study the behavior of higher chiral-oder terms and to examine a possible 
connection between these higher order terms and the heavy-meson exchange contributions, 
it seems of great importance to perform a v = 2 calculation. Such a calculation has 
recently been done by Dmitrasinovic, Myhrer, Sato and myself (DKMS) p9[ . I describe 
here briefly the highlights of our results. 

As an initial attempt we may concentrate on the effective transition operators 
themselves instead of the full distorted-wave transition amplitude T. Furthermore, we 
limit ourselves to the threshold kinematics, which means that, in Fig. 1(b), q = (m n ,0) 
and k = (m n /2, k) with k 2 = —m^m^. Enumerating u—2 irreducible diagrams that give 
rise to transition operator T^ =2 ) for NN — > NNir, we find 20 topologically distinct types 
of diagrams (nineteen of them are new). For a particular case of the pp — > ppn° reaction 
near threshold, the isospin selection rules and the s-wave character of the outgoing pion 
reduce this number from 20 to 7 (six of them are new). Some representative diagrams are 
depicted in Fig. 2. 





(a) (b) (c) (d) 

Fig.2: Selected u — 2 diagrams; the thick (thin) lines represent nucleons (pions). 

Let T^ u=2 ^ (n = 1,2 ... ,7) represent the transition operator coming from the n-th type 
of diagrams. The importance of each operator may be measured in terms of the ratio 
1Z n = T( v=2 > /T( u=l \ where T^ u=1 ^ is the leading-order two-body transition operator. 

Of the seven types of diagrams, some can be interpreted as vertex corrections to the 
lower order operators, T\ v= ~ x > and T^ v=l >\ Fig.2 (a) gives an example. The other types are 
rather loosely called "two-pion" diagrams; three examples are shown in Figs. 2(b),(c),(d). 



We have found [29 that TZ^s corresponding to the vertex-correction-type diagrams are 
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small (0.1 ~ 0.2) in conformity with the general expectation of xPT. On the other hand, 
some of the "two-pion" diagrams turn out to give very large contributions. Especially, 
TZJs belonging to the types illustrated in Fig. 2(b), Fig. 2(c) and Fig. 2(d) are individually 
very large. For the first two, \7Z n \ = 3 ~ 7 (depending on the input low-energy parame- 
ters), while \TZ n \ = 7 ~ 10 for the pion-pion rescattering diagram, Fig.2(d). This feature 
is consistent with the expectation that the pp — ► ppir reaction is sensitive to heavy- 
meson exchanges. Thus the phenomenologically important a-meson contributions p3| 
seem to have discernible "representatives" among the NNLO xPT diagrams. The large 
contributions of the "two-pion" diagrams do not necessarily constitute evidence for the 
non- convergence of xPT expansion. Since these diagrams appear only at NNLO or higher, 
the convergence can be tested only by calculating corrections to the NNLO diagrams. On 
the historical note, I should mention that Gedalin et al. 2£| considered some of the NNLO 
diagrams and pointed out that they could be large. According to our calculation, some 
important NNLO diagrams are missing in |28) and a loop integral expression in |28 needs 
to be corrected. 

DKMS use the "standard" chiral counting rule of Weinberg 0, whereas Cohen et al. 
EBf have emphasized that for the NN — > NNtc reaction, which involves significant energy- 
momentum transfers, one should modify the chiral counting rule. In this new counting 
scheme, the expansion parameter is not any longer m^/m^ but y^m^m^. It is a future 
task to perform a calculation similar to that of DKMS using the modified counting rule. 
Furthermore, to obtain transition amplitudes that can be directly compared with the 
experimental cross sections, it is imperative to carry out a distorted-wave calculation. We 
hope to be able to report our investigation along this line in the near future. 

In the first part of my talk I have surveyed the recent progress in \PT treatments of 
low energy-momentum transfer observables in the two-nucleon systems. I have mentioned 
that we are entering the first stage of ab initio calculations based on HB^PT. It will 
be extremely nice if we can perform similar ab initio calculations for e.g. NN — > NNw, 
which involve higher energy-momentum transfers. This is a difficult but urgent challenge. 
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